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INTRODUCTION 


* In 1934 , Gotz and Dobson made a fundamental study on the 
problems associated with inversion of radiance measurements 
for atmospheric ozone profiles. 


* First method for vertical atmoscr; 
from satellite measurements prese 
Wentworth . 
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* The current National Space Science Data Center algorithm 
for BUV (Nimbus-4) and 3BUV ; Nimbus-"; archives is base* 
primarily on the work of C.D. Rodgers. 


* This method formulated m terms error covariance 
matrices, is associate! with both direct .u- lourements . 
apriori first guess profiles, and a weighting function 


matrix . 


STATEMENT OF THE PROBLEM 


Profile Retrieval 
Power Spectrum 


1. Profile retrieval: 
Transfer Equations. 




Achieved by Inversion of Radiative 

c/v - 


; Fredholm Integral/lst kind) 


K = Kerne 1 1 PI ar.rk Radiation Law/Maxwell Equations) 
x = unknown profile 

y = given data set of radiance values 

* To numerically solve: K(x) assumed independent of x. 
Then the equation is written in a linearized form: 

y = Ax f >£ 


A = weighting function matrix 

— HO 1 36 


* Standard Approach 
-- Ignore 
— Find A 1 
— Solve x = A 'y 


* Prob 1 ems : 

1. The A matrix is near singular 

2 Number of Unknowns in the profile exceeds the 
number of data points (The Ill-Posed Problem) Therefore, 
the solution may not be unique. 

3. Even if a unique solution exists. ^ may cause 
the solution to be ill-conditioned. 


maximum entropy solution 

* Since the number of unknowns exceeds the number uf “u^a 
points, probability theory is needed. 

* Introduce Maximum Entropy forma 1 ism which induces an 
unknown probability distribution from partial data. 

* Maximize the information measure subject to the 
following constraints: 

Sum of the probabilities is 1. 

The data, given as averages, is written in the 
form of expectation values. 

* Solution: An exponential probability distribution. 

p t = exp [ - i A » . ] / <- ) 

where = Lagrange multipliers, and 

= Partition function. 

2. Power spectral estimation for a time series of TOMS 
data . 

In 1967 Burg introduced his Maximum Entropy Method (MEM) 
for power spectral estimation 

(Pn+ U/fN 


P(0 
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2 H + £ a N « 1-2* | 2 

n=l 

- P N is the prediction filter error power 

- Snm are the filter coefficients 

- frg is the Nyquist/sampl ing frequency 

The coefficients ere obtained from the Yule-Walker 
equat ions : 
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- are the time -core 1 at l ons functions or the 

_ Length of the filter. N. has to be determined. 


RESULTS 


The cases depicted in the Figures are preliminary 
and be ing^re fined for publication m two separate journals 
of the joint American and European physical - 
ngures la and lb show the case of the o:one f 
retrieval as done by then National opace oCi *v ce ,^ a Th 
Center end the Maximum Entropy Met ^ «sp a * ™ 

latter shows a clear depletion _o. --one or ^ 

STM MSSUMpUnS: the Maximum Entropy 

Shows a more realistic ozone depletion P-c-ie. 

p ig „ v - e q n a 3b and 3c involve the time-series and 

r igu* e- - ^ city of Hyderabad. 

Maximum Entropy ^ 1 ; * The randomness and 
India for the ac_i ~ T V, M r' t :rne series is clearly 

periodicities eviue ~ or r ant . the three expected 

confirmed in Figure ^ -me _ le _ the mid-year 

abnormalities, name and the annual cycle of 

peax du T -v th .... - ar p present. This is strong 

the motion of the en. a it-* - dl - i- rc '=''“ Th , ther 

evidence for the validity of this approach. The other 

cvcles are being studied with perhaps a correlation or tne 
210 day cycle with the half period of the chandler wcbiwi<r 
as a possibility. The 80 day sample was performed as a 
check of both the stationary of tne series and tne 
authenticity of the 365 spectrum cycles (rig ~»c ) . This 
confirms the origin of the physical processes governing the 
spectral periodicities as arising very eariy within the 
solar year! 

Figures 4a and 4b are the Id MO time series and full 
Maximum Entropy spectrum for Hyderabad for the less active 
1980 year. A slight shift in both the dynamic range and 
frequency again cleany show that the ^er^es are not quite 
stationary from year to year, but the major cycles remain. 
In all cases, a comparative FFT ope ct rum cua.d resolve no 
more than the mid-year peak and was therefore abandoned. 

The remaining cases show both the time series and 
resulting full and partial sample spectra New York City 
(41* lat ) and SYOWA < Japanese ) antarctic station ( /u iaii. 
The New York data shows a much -larger random spread than 
Hyderabad and this resulted in a larger number of multiple 
frequencies in the spectrum. A greater number of major 
cycles are present as are some of those found in Hyderabad 
data including the speculative half period Chandler wobble. 

The SYOWA spectrum was calculated r or the cases of 
134 days (up to the Sun gap), and 365 days ( inc luding the 
gap) and again, a short sample of 80 days. This was done 
in part to study effects such as the quasi- oi 1 mear- 
oscillation (QBO) . More work needs to be done at this 
point to estimate correctly the secular trends in tota 
ozone . 
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